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ABSTRACT 

We present new FUSE spectroscopy of the dwarf starburst galaxy NGC 625. 
These observations probe multiple phases of the interstellar medium, including 
the coronal, ionized, neutral and molecular gas. This nearby (D = 3.9 ± 0.2 Mpc) 
system shows a clear detection of outflowing coronal gas as traced by O VI A 1032 
A absorption. The centroid of the O VI profile is blueshifted with respect to 
the galaxy systemic velocity by ~ 30 kms~^, suggesting a low- velocity outflow. 
The implied O VI velocity extent is found to be 100 ± 20 kms~^, which is 
fully consistent with the detected H I outflow velocity found in radio synthesis 
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observations. We detect multiple lines of diffuse H2 absorption from the ISM of 
NGC 625; this is one of only a few extragalactic systems with FUSE detections 
of H2 lines in the Lyman and Werner bands. We find a potential abundance 
offset between the neutral and nebular gas that exceeds the errors on the derived 
column densities. Since such an offset has been found in multiple dwarf galaxies, 
we discuss the implications of a lower-metallicity halo surrounding the central star 
forming regions of dwarf galaxies. The apparent offset may be due to saturation 
of the observed O I line, and higher S/N observations are required to resolve this 
issue. 

Subject headings: galaxies: starburst — galaxies: dwarf — galaxies: abundances 
— galaxies: kinematics and dynamics — galaxies: individual (NGC 625) 



1. Introduction 

Local starbursting dwarf galaxies offer a unique opportunity to study the various phases 
of the interstellar medium (ISM) under extreme conditions and in unprecedented detail. 
Furthermore, the physical conditions in these targets are likely very similar to those found 
in galaxies at high redshifts, making them crucial benchmarks against which observations 
of such targets must be compared. With these observations, we capitalize on the high- 
resolution capabilities of the Far Ultraviolet Spectroscopic Explorer {FUSE; see Moos etal. 
2000) to observe the interaction of the multiple phases of the ISM in NGC 625, a nearby (D 
= 3.9±0.2 Mpc; Cannon etal. 2003), low mass (Mm = 1.1x10^ M©; Cannon etal. 2004), 
low metallicity (12 + log(0/H) = 8.14 ±0.02; Skillman etal. 2003b) dwarf irregular galaxy 
in the Sculptor Group that is currently undergoing a massive star formation episode (see 
Table 1 for basic galaxy properties). 

Of paramount importance in the evolution of dwarf galaxies is the role played by out- 
flows powered by the energy inputs of supernovae and stellar winds. These outflows may 
drive significant fractions of the gas from galaxies with shallow potential wells. Recent simu- 
lations suggest that the return of metals into the intergalactic medium (IGM) finds its most 
effective avenue in these low mass galaxies (Mac Low & Fcrrara 1999). Understanding the 
nature of these outflows is therefore a fundamental astrophysical problem with far-reaching 
imphcations. These FUSE spectra isolate the powerful diagnostic absorption hues of O VI 



^ Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE 
is operated for NASA by the Johns Hopkins University under NASA contract NAS 5-32985. 
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AA1032, 1038 A which probe the temperature and kinematic structure of the coronal gas 
in a galaxy (see, e.g., Hcckman etal. 2001). The O VI ion traces hot ionized gas over the 
temperature range T ~ 10^ - 10^ K, with maximum sensitivity to gas in collisional ionization 
equilibrium at temperatures T ~ (2 - 3)xl0^ K (Sutherland & Dopita 1993). Thus, knowing 
the temperature of the emitting material, we can (with some assumptions) directly estimate 
electron density and pressure in the coronal gas. We can then assess the importance and 
timescale of radiative cooling in outflows and in the starburst phenomenon. Indeed, in a 
case such as NGC 625 where the outflow is of comparatively low velocity (see, e.g., Heckman 
etal. 2001), such effects may be dominant. These observations thus provide an opportunity 
to investigate many properties of outflows in dwarf galaxies. 

Multiple independent lines of evidence suggest that the recent star formation in NGC 625 
has been violent and that the coronal gas content is substantial. An extended soft x-ray 
component was detected in RO S AT ima,gmg (Bomans & Grant 1998) and has been verified 
in newly-obtained Chandra imaging (Cannon etal., in preparation). To attain hot gas at 
large distances from the central starburst requires the presence of an outflow in the recent 
past. This outflow is also seen in H I synthesis imaging (Cannon et al. 2004), making NGC 625 
one of only a few dwarf galaxies with a detected neutral gas outflow. This outflow appears 
to be a result of the extended (both spatially and temporally) star formation event that 
the galaxy has undergone over the last ^ 50 Myr (see Cannon etal. 2003 for details). This 
extended burst is unexpected, given the brevity of star formation implied by the presence of 
the A 4686 A spectroscopic Wolf-Rayet (W-R) feature (^ 6 Myr; see Conti 1991; Schaerer 
etal. 1999; the spectrum is presented in Skillman etal. 2003b). 

Theoretical and observational evidence suggests that star formation behavior should be 
a function of metallicity (see, e.g., Maloney & Black 1988; Taylor etal. 1998). NGC 625 
has a metallicity significantly below the solar value, [0/H] = —0.47 ±0.06 (Skillman etal. 
2003b), which provides an opportunity to study the interplay between powerful starburst 
activity and molecular gas content at a metallicity similar to that of the Small Magellanic 
Cloud. Single-dish CO spectra reveal a large CO cloud coincident with the main dust 
concentration (Cote, Braine, & Cannon, in preparation), but offset from the current massive 
star formation regions. This position also hosts a moderately-extincted, low-mass stellar 
cluster that produces a thermal radio continuum peak (Cannon & Skillman 2004). 

In light of these arguments, we also analyze the diffuse H2 content of this relatively 
low- metallicity galaxy. The Lyman- Werner bands of II2 arc expected to be very sensitive 
to relatively cool diffuse H2 gas (T ~ 100 - 1000 K) over a wide range of column densities 
(N(H2) ~ 10^^ - 10^^ cm~^). In observations of galaxies or large star-formation regions, these 
fines probe the diffuse H2 content along many sightfines to ultraviolet (UV)-luminous sources. 
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Evidence suggests that diffuse H2 clouds are indeed prevalent at low metallicities; Tumlinson 
et al. (2002) detect diffuse H2 along 92% of sightlines toward UV sources in the SMC, which 
has very nearly the same metallicity (Dufour 1984; Garnett 1999) as NGC625. However, 
integrated FUSE spectra of metal-poor starbursts suggest that diffuse H2 is more difficult to 
detect in distant targets (Vidal-Madjar etal. 2000; Heckman et al. 2001; Thuan etal. 2002; 
Aloisi et al. 2003; Lebouteiller et al. 2004; Lecavelier des Etangs et al. 2004; Hoopes et al. 
2004; see further discussion in § 4.2). With these data we add NGC625 to the small but 
growing sample of dwarf starburst systems with FUSE observations of diffuse H2 gas; we 
detect multiple H2 lines and discuss the implied properties of the molecular gas. 

Finally, since the spectral region probed by FUSE is rich in neutral and ionized gas 
absorption lines, we also address the column densities and abundances of observable species 
in the ISM of NGC625. This point is especially important since evidence is growing for 
a bimodal abundance distribution in the ISM of dwarf galaxies. The nebular regions (i.e., 
near the star formation regions and accessible to abundance studies via optical and near- 
infrared emission line spectroscopy) appear to have elevated abundances compared to the 
neutral interstellar gas through which FUSE sightlines usually pass. The sample of dwarf 
systems with such an analysis is small, but the offsets appear to be pronounced (> 0.5 dex in 
oxygen abundance). With these arguments in mind, we analyze the neutral gas abundances 
in NGC 625 and compare these to values obtained from nebular spectroscopy obtained by 
Skillman etal. (2003b). 

Understanding the kinematics and behavior of the outflowing coronal gas, the molecular 
gas content, and the neutral-gas abundances are the major goals of this work. With these 
data, we can characterize the interaction of the outflowing gas with the neutral components as 
revealed in our H I imaging and from other ionization species within the FUSE spectrum. In 
addition, they allow us to estimate the metallicity of the neutral gas and compare these values 
with those found for emission line gas in the nebular regions. We discuss the evolution of 
dwarf starburst galaxies using these FUSE data and published information in the literature. 

2. Observations and Data Reduction 

FUSE spectroscopy of NGC 625 using the 4" x 20" medium resolution aperture (MDRS) 
was obtained on 2003, November 8-9, for program D040. No roll angle constraints were placed 
on the observation, since the aperture width is comparable to the size of the high-equivalent 
width Ha emission (and hence the most UV-luminous sightlines). Our observations were 
obtained with an average slit angle 112° east of north; the approximate aperture placement 
is shown superposed on Hubble Space Telescope (HST) WFPC2 V-band and (continuum 
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subtracted) Ho; images in Figure 1. The total integration time on-source was 58.1 ksec, 
with ~ 76% of this occurring during orbital night, which decreases the contamination of the 
spectra by terrestrial N I and O I airglow lines. Data from both orbital day and night were 
used in spectral regions not affected by airglow lines; otherwise, only data obtained during 
orbital night were used in our analysis of the NGC 625 absorption features. 

This observation produced 23 raw time-tagged exposures which were re-processed using 
a recent version of the FUSE calibration software (CALFUSE v2.4) available from Johns 
Hopkins University^. This pipeline reduction removes mirror, grating, and spacecraft mo- 
tions, and then corrects for astigmatism and Doppler motions. Wavelengths are assigned to 
each photon, and these events are screened for data quality (e.g., event bursts, spacecraft 
jitter, etc.). The data are then flux calibrated and the bad pixel maps are corrected for 
spacecraft motions. Finally, spectral extractions are performed for each detector segment for 
each channel using either the orbital night-only data or the orbital day and night data. Each 
individual exposure was processed by the pipeline and the resulting spectra were combined 
to produce the final spectrum for each channel and data segment. This results in a total of 
16 extracted spectra. 

FUSE uses four optical channels to produce 8 detector segments that allow spectral 
extractions between ~ 900 and 1200 A. The two SiC channels are optimized for shorter 
wavelengths (905 to 1100 A), while the LiF channels are optimized for longer wavelengths 
(980 to 1187 A). The LiF channels are more sensitive than the SiC channels, so these provide 
our highest signal-to-noise (S/N) data. The effective area of the telescope is also maximized 
at ~ 1032 A (i.e., very near the important O VI absorption features), so we concentrate our 
analysis on the LiFl channel, with LiF2 (lower S/N) serving to verify detections. Data in the 
SiC channels are generally of low S/N (caused by a combination of lower effective area and 
channel misalignment between the LiF and SiC channels), and we do not use them in this 
analysis. We present in Figure 2 the LiFlA and LiFlB spectra of NGC 625, binned to 0.05 
A resolution for clarity. Note the richness of the spectrum, with numerous strong absorption 
lines detected. We do not combine any of the overlapping detector segments, due to the 
changing spectral resolution and sensitivity of the detectors as a function of wavelength. 
The final velocity resolution is ~ 30 kms~^, with S/N ratios of ~ 10 and 5 per resolution 
element at 1032 A for the LiFl and LiF2 channels, respectively. All velocities quoted are in 
the heliocentric reference frame. 



^See http://fuse.pha.jhu.edu/analysis/calfuse.html 
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3. Analysis 

These FUSE data probe four different ISM phases: diffuse molecular gas, neutral gas, 
warm photoionized gas, and hot coronal gas. We discuss the kinematics of each of these 
phases as derived from Gaussian fitting to the absorption profiles in the following subsec- 
tion and summarize important line properties in Table 2. We then discuss derived column 
densities of the gas in § 3.2 and the inferred elemental abundances in § 3.3. 

3.1. Gas Kinematics 

There are two main absorbing components seen in these data. First, Milky Way gas 
appears near zero velocity. Second, NGC625 absorptions are centered around + 400 kms~^. 
In addition, there may also exist a halo intermediate- velocity cloud as seen in the red wings 
of some of the stronger Milky Way absorption lines (see below) . In this work we concentrate 
on the lines arising within NGC625. 

Beginning with the molecular phase, we detect numerous low-level H2 absorption lines 
in NGC625. The Lyman and Werner bands of H2 produce hundreds of absorption hues 
in the spectral region under study. The lowest rotational levels (J = 0, 1, and 2) are 
expected to be most highly-populated by the moderate diffuse molecular hydrogen columns 
and temperatures typically found in external galaxies (N(H2) ~ 10^^ - 10^^ cm~^; T ~ 100 
- 1000 K). Of these low-level lines, six are detected in unblended, high S/N regions of the 
final spectra. We do not detect any lines in the higher J levels. The velocity centroid of 
the detected lines (including R(0)(2-0) A 1077.1399 A, R(l)(4-0) A 1049.960 A, R(2)(5-0) 
A 1038.689 A, P(l)(5-0) 1038.157 A, R(2)(0-0) A 1009.024 A, and R(l)(8-0) A 1002.449 A) is 
found to be vh2 = 401 ±10 kms~^ (all errors quoted in this work are la unless otherwise 
noted). This is in general agreement with, though slightly lower than, the H I systemic 
velocity derived by Cannon etal. (2004), W^ys = 413 ±5 kms~^ We note that our FUSE slit 
placement covers H I gas in emission at velocities from ~ 405 - 415 kms~^. Since velocity 
offsets of up to 20 kms~^ are common in FUSE spectra, we interpret the velocity centroid 
of the molecular gas to be coincident with the H I and stellar populations. We measure 
blueshifts of other species with respect to the velocity centroid for H2 absorption, vhj — 
401 ±5 kms-^ 

We detect absorption lines of N, O, Si, P, Ar and Fe that arise from the neutral gas 
phase (i.e., ionization potentials greater than or near that of H). The velocity centroids of 
the strongest lines in clean spectral regions (including N I A 1134.980 A, N I A 1134.415 A, 
O I A 1039.230 A, Si II A 1020.699 A, P II A 1152.818 A, Ar I A 1048.220 A, Ar I A 1066.660 A, 
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Fe II A 1096.877 A, and Fe II A 1063.176 A) yield an average neutral gas velocity of Vneutrai = 
392 ± 10 km s~^. This can be compared with the molecular gas velocity derived above, Vsys = 
401 ± 10 kms~^. These values are equal within the errors; the small offset between the two 
may be caused by the neutral gas being split into two kinematic components (separated by 
10-20 kms~^). To test for this, we compared fits for single and double Gaussian components 
to the strongest line profiles (including O I A 1039.230 A, Si II A 1020.699 A, and Fe II 
A 1063.176 A), but found no statistically significant difference (in per degree of freedom) 
between them. The weakness of the putative second component introduces a negligible 
contribution to the error budgets for derived column densities and abundances (see below). 

There are two ionized gas absorption lines detected in clean regions of our spectra, 
including the A 1037.018 A line of C 11* and the A 1012.495 A line of S ill. These absorptions 
give an average ionized gas velocity of vionized — 377 ±9 kms~^; i.e., offset with respect to 
the molecular gas by —24 kms~^. Finally, coronal gas is detected in O VI A 1031.926 A 
absorption at a velocity of 370 kms~^. This is blueshifted from the molecular gas velocity 
by —31 kms~^, from the neutral gas by —22 kms~^, and from the ionized gas by —7 kms~^ 
(though equal within errors). The weaker line of the doublet at A 1037.617 A falls in a 
complicated spectral region and is not used in the present analysis; see Figure 3 for a closer 
view of the spectral region around O VI. We present in Figures 4 and 5 the normalized 
hue profiles of the detected absorption lines from the ISM of NGC 625 (absorption features 
were normalized by using clean spectral regions surrounding the lines; our error budget 
incorporates errors as a result of the continuum placements, which are usually comparable 
to the statistical noise) . 

This scenario is consistent with a low-velocity outflow, where the coronal gas has the 
largest velocity offset from the starburst region (which presumably is coincident with the 
velocity of the molecular and neutral gas absorption features). However, it is important to 
note that although the velocity centroids of the different types of gas may differ, there is 

considerable overlap in the velocity extents of many of the profiles (examine Figures 4 and 
5). We interpret these changing velocity centroids and line breadths as results of the outflow, 
but it is clear that gas of all ionization levels co-exists throughout the galaxy. This velocity 
gradient is discussed further in § 4.1. 



3.2. Column Densities 

For each unblended spectral line shown in Table 2, the absorption proflles were converted 
into apparent optical depths, and then these proflles were integrated over velocity (see further 
discussion in Savage & Sembach 1991). The essence of this technique is to measure the depth 
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of normalized absorption lines as a function of velocity, and then to use atomic physics to 
infer the column density that best reproduces such a line profile. The column density and 
normalized line depth are related by: 



where N is the column density in atoms cm~^, mg is the mass of the electron, c is the speed 
of hght, e^ is the standard unit of charge, / is the oscillator strength of the transition, A is 

the rest wavelength of the transition, and In^j^ is the apparent optical depth as a function 
of velocity (the natural logarithm of the ratio of intensity in the continuum versus intensity 
in the line, as a function of velocity). This method is preferred over other common practices 
(e.g., curve of growth analysis) because it docs not presume a functional form for the line 
under analysis. In general, the results from standard curve of growth techniques and the 
apparent optical depth method agree very well, and we adopt the apparent optical depth 
method values in this work. Note that the kinematic characteristics were derived separately 
using Gaussian fitting (see above). 

For the diffuse molecular hydrogen absorption detected, we derive weighted mean col- 
umn densities of log(N(H2)) = 14.92±0.19 (J = 0), 14.91±0.09 (J = l), and 14.47±0.13 
(J = 2). Following Tumhnson etal. (2002), we can infer a range of rotational temperatures, 
Toi, of ~ 70-90 K. At high densities where collisions dominate the level populations, this 
value is indicative of the kinetic temperature of the gas. Such a progression is expected for 
relatively low-excitation H2, with the bulk of the diffuse gas in the lowest available rotational 
levels, and the column densities rapidly decreasing toward higher rotational levels. Deriva- 
tion of the H I absorbing column is discussed in § 3.3; combining this estimate with the 
measured molecular columns implies a diffuse molecular fraction ~ 2xl0~^. It should be 
noted that this value samples the most UV-bright sightlines, and hence is heavily weighted 
toward material on the line of sight to the major starburst region. This is not coincident 
with the highest H I column density seen in emission (which is offset with respect to the 
starburst region; see Cannon et al. 2004) nor with a large detected molecular cloud (from 
CO observations; see Cote, Braine, & Cannon, in preparation). In general this highlights a 
shortcoming of the Lyman and Werner bands in probing cool molecular gas, as these lines 
are sensitive to diffuse gas (i.e., sightlines with low extinctions) but are less sensitive to the 
more common discrete molecular clouds and star formation complexes where the bulk of 
the cool molecular material is expected to reside (since these regions will be more heavily 
extincted and hence have less background UV continuum emission). This point is discussed 
in more detail in § 4.2. 




(1) 
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For other neutral and ionized gas absorption lines present in our spectra, we derive the 
following column densities: log(N) = 14.63 (N l), 15.80 (O l), 15.37 (Si ll), 13.54 (P ll), 
15.04 (S III), 13.98 (Ar l), and 14.73 (Fe ll); typical errors are 0.1 - 0.2 dex (see Table 2). 
We note that the O I A 1039.230 A line may be saturated; since the SiC channels are of low 
S/N we cannot probe the amount of saturation empirically with other O lines in these data 
alone. We discuss tests for saturation and their implications in § 3.3. For the purposes of 
this paper, we adopt the empirical O I column (15.80 ±0.12) as a lower limit; we note that 
this implies an upper limit for our implied neutral gas [N/0] ratio (see next section). For the 
coronal gas as traced by O VI A 1031.926 A absorption, we find log(N) = 14.32 ±0.08. We 
use the neutral species to derive gas-phase abundances in the next subsection; the coronal 
gas is discussed in more detail in § 4.1. 

3.3. Relative Gas-Phase Abundances 

3.3.1. The Intrinsic Hi Column 

An estimate of the H I column density in NGC 625 from these data alone is hindered 
by the low S/N of the spectra in the SiC channels; this precludes us from using the higher- 
order Lyman lines, and we are left with only Ly/3 to constrain the H I column. However, 
the velocity separation of NGC 625 and the Milky Way and the low measured Galactic 
column along this high-latitude sighthne (—73.1°; we find a Galactic column log(N(H l)) = 
20.3 ± 0.2) are sufficient to clearly separate Galactic and intrinsic absorption profiles. This 
allows us to analyze both the red and blue wings of the NGC 625 H I absorption feature, 
better constraining our fit to the H I column density. We present in Figure 6 a closer view of 
the region surrounding Ly /3, overlaid with fits to both the Galactic and intrinsic H I columns. 

A pronounced uncertainty remains in our measured H I column due to the presence of 
strong stellar O VI P-Cygni profiles in the spectrum redward of Ly/3 (see Figure 6). These 
features result in an uncertainty in the adopted continuum in the region around Ly/5. We 
adopt the mean of a "low" and "high" continuum placement, as shown in Figure 6. This leads 
us to estimate the intrinsic H I column in NGC 625 at log(N(H l)) = 20.5 ± 0.3. This is at least 
half a dex lower than the H I 21 cm emission column density in this region; comparing the 
FUSE aperture placement (see Figure 1) with the ~ 15" resolution H I column density map 
presented by Cannon etal. (2004), we find that the FUSE aperture is completely enclosed 
within the 10^^ cm~^ contour. We believe that this estimate and associated uncertainty for 
the sight line probed are robust. Using the 21 cm value would result in a stronger Ly /3 profile 
than the other models shown in Figure 6 (see the dot-dash fine corresponding to log(N(H l)) 
= 21). The difference between our absorption estimate and the 21cm value could easily be 
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due to the much different sohd angles probed by the two measurements (a size ratio of 2.2:1 
exists between the H I emission beam and the sohd angle of the FUSE MDRS aperture). We 
discuss the implications of this offset in H I absorbing column in more detail in § 4. 

3.3.2. Neutral Gas Column Densities and Implied Abundances 

Using these column density estimates and assuming that the observed absorption lines 
arise from the main elemental species in the neutral ISM, we find the following abundances 
and elemental ratios: [N/H] = —1.80 ±0.34; [0/H] > —1.36 ±0.33 (wc stress again that the 
O I A 1039.230 A line may be saturated, and hence the lower limit is quoted; see further 
discussion in § 3.3.3); [Ar/H] = -1.08 ±0.33; [Si/H] = -0.67 ±0.34; [Fc/H] = -1.22 ±0.34; 
log(N/0) < -1.17±0.17 ([N/0] > -0.44±0.21). These values use the updated solar 
abundances of Holweger (2001; N, Si, Fe) and Asplund et al. (2004; O), or the older standard 
values of Anders & Grevesse (1989; Ar). These abundances are summarized in Table 3. Note 
that a downward shift of the assumed H I column leads to increased abundances, and that 
an increase in N(H l) leads to a decrease in inferred metal abundance (specifically, if the H I 
column is a factor of 2 larger, as might be inferred from the 21 cm emission in this direction 
[see § 3.3.1], the derived abundances will be a factor of 2 lower than the quoted values). 
Also, the N/0 abundance ratio assumes that the neutral gas components are co-spatial (i.e., 
that the ionization correction factor is negligible). 

Keeping in mind the relatively large errors involved with these absorption line abundance 
measurements, and the points about potential line saturation discussed in § 3.3.3, these 
neutral gas values can be compared with the nebular abundances presented by Skillman 
etal. (2003b). Therein, for the major star formation region which dominates these FUSE 
sightlines, the calculated abundance ratios were found to be [0/H] = —0.47 ±0.06 and 
log(N/0) = —1.33 ±0.02. Our neutral gas oxygen abundance (taken at face value) is found 
to be lower by ~ 0.9 dex, while the relative abundances between N and O are very similar. 
The offset between the neutral and nebular gas abundances may indicate that the outer 
regions of dwarf galaxies are less enriched in heavy elements compared to the inner, star- 
forming regions. We discuss this point further and compare to results found for other dwarf 
starbursts in the literature in § 4.3. 

It is interesting to note that Skillman etal. (2003b) found an elevated N/0 ratio in the 
three highest surface brightness H II regions in NGC 625 when compared to other star-forming 
dwarf galaxies. This vahic was found to be comparable to those for blue compact dwarfs, 
and was postulated to potentially arise as an effect of a long quiescent period prior to the 
current (temporally and spatially extended) star formation episode where the N abundance 
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may be elevated by the evolution of intermediate mass stars. Our neutral gas N/0 ratio 
(which is equal to the nebular value, within errors) is also found to be high compared to 
the neutral N/0 values derived for other dwarf starbursts observed with FUSE (note that 
if the O I A 1039.230 A line is saturated, the N/0 value decreases, moving the NGC625 
N/0 ratio closer to the values found for other systems). While interpretation of the elevated 
N/0 ratio in the neutral gas of starbursting dwarfs will require a larger statistical sample, 
it appears that in NGC625 the nucleosynthetic products of intermediate mass stars have 
been enriching the ISM for a long period, and that these products have been mixed into the 
neutral gas at large separations from the starburst regions (see also § 4.3). 



3. 3. 3. Effects of Potential Line Saturation 

As mentioned previously, the O I A 1039.230 A line may be saturated, potentially severely 
underestimating the true O I column density. As shown in Pettini & Lipman (1995), strongly 
saturated O l absorption lines can imply abundance ranges of up to ~ 1000 without differ- 
ences in the goodness of the profile fit. Examining the profile of the O I line in Figure 5, it is 
clear that even this moderate-strength transition [log(/-A) = 0.974] is close to saturation in 
the central regions of the line. Furthermore, since these FUSE daia, probe multiple sightlines 
(see Figure 1), the line may be more heavily saturated than the profile shows, if a luminous 
continuum source is not shielded by the foreground O I region. For these reasons, we more 
carefully explore the tests for, and potential effects of, O I line saturation. 

We can probe the amount of line saturation by subtracting the local rms value from the 
line profile and re-deriving the column density via the apparent optical depth method. This 
exercise shows that the center of the line is within 1 cr of zero fiux; this causes the inferred 
column density to grow quickly (see Equation 1 and discussion in § 3.2). This suggests that 
the O I fine is at least partially saturated. 

The degree of this saturation can be tested by using tracer species to infer the column 
of neutral oxygen. First, Si II can be used as a proxy by assuming that the 0/Si ratio is 
the same as the solar ratio (see Lu etal. 1998 for details). While uncertain, this method 
suggests that O l saturation may cause the O I column to be underestimated by as much as 
0.7 dex. Note, however, that Si can be depleted onto dust grains, and thus this technique is 
not without error. Second, P ll (which is less depleted onto dust grains) can be used as a 
tracer of the neutral oxygen abundance. Following Lebouteiller etal. (2004), if [P/0] = is 
assumed for the neutral gas, then the implied column of O I is nearly a full dex above that 
derived from the actual O I line profile. Both of these proxy methods imply that the O I 
column has indeed been severely underestimated. 
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Ideally, one would use other oxygen lines in the SiC channels to independently test for 
line saturation. As mentioned previously, these data lack the S/N in the SiC channels to 
perform such a test. Thus, in lieu of a large sample of 0/Si and P/0 studies of the neutral 
ISM of dwarf galaxies, for the remainder of this paper we explore the results of a column 
of O I as derived above, log(N) = 15.80 ±0.12, which implies an abundance offset when 
compared to the nebular regions. We discuss the implications of this scenario in § 4.3. We 
stress, however, that this O I column density is a lower limit and that the effects of hne 
saturation may indeed be pronounced. 

4. Discussion 

The most important results in this paper arc as follows: 1) The discovery of outflowing 
coronal gas from a relatively low-luminosity dwarf starburst galaxy; 2) The discovery of 
diffuse H2 absorption in a relatively distant, low-metallicity galaxy; 3) The discovery of a 
possible abundance offset between the nebular and neutral gas phases in a star-forming dwarf 
galaxy. In this section we discuss each of these points in more detail and compare the values 
found for NGC 625 with those found for other dwarf starbursts in the hterature. 



4.1. Outflowing Coronal Gas 

The coronal gas content of low-mass systems that are undergoing outflow episodes is an 
important component of models of superbubble evolution. If the coronal gas density is high 
and radiative cooling is efficient, then outflow energies can be expended prior to the outflow 
breaking out of the disk and venting hot gas, metals and energy into the surrounding IGM. If, 
on the other hand, the coronal gas density is low and the radiative cooling is inefficient, then 
the coronal gas does not radiate sufficient energy to slow or stall galactic-scale outflows from 
low-mass galaxies. While the evolution of outflows is complex and depends on various other 
factors, coronal gas cooling is one important effect that has only recently become accessible 
observationally. 

In principal, by observing the absorption and emission profiles of outflowing O VI gas 
from a galaxy, one can obtain an estimate of the cooling rate of the gas and compare this to 
the estimated energy input into the ISM by the evolution of the massive stars associated with 
recent star formation (see Heckman etal. 2001). However, such an estimate of the cooling 
rate in NGC 625 is difficult for two reasons. First, it is not straightforward to extrapolate to 
an unambiguous outflow cavity size estimate. Second, the low outflow velocity of the O VI 
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gas with respect to the molecular gas (and presumably the stellar population) in NGC 625 

(~ 30 kms~^) does not cleanly separate the absorption profile from the expected location of 
O VI emission (if we see emission only from the redshifted, back side of a symmetric outfiow 
along the line of sight). These factors complicate estimates of the mass and density of the 
O VI gas, both of which are needed to constrain the cooling rate (see ShuU & Slavin 1994). 

FUSE observations of O VI emission from other starbursting galaxies have shown that 
typical cooling rates in coronal gas are small. Heckman et al. (2001) and Hoopes et al. (2003) 
estimate that coronal gas cooling is roughly comparable to the energy radiated away in 
soft x-ray emission in the outfiows in NGC 1705 and M82. Since the soft x-ray emission 
radiates only a small fraction of the input energy of the wind (Strickland & Stevens 2000), 
the coohng rate via coronal gas appears to be small. A similar column of O VI gas has 
been found in NGC 625 (log(N(0 VI)) = 14.32 ±0.08) and NGC 1705 (log(N(0 VI)) = 
14.26 ±0.08); because the outflow cavity is hkely smaller in NGC 625 than in NGC 1705 
(implying a smaller mass of cooling plasma), it is likely that radiative cooling via coronal 
gas is not a significant portion of the energy budget of NGC 625. 

Comparing to the literature, it appears that pronounced O VI absorption is ubiquitous 
in actively star-forming systems with non-zero inclinations. Extensive observations through 
the Milky Way halo have shown an average O VI absorption column density of 14.38 (Savage 
etal. 2003). Similarly, Howk etal. (2002) and Hoopes etal. (2002) have shown that O VI ab- 
sorption is present along tens of sightlines into the Large and Small Magellanic Clouds, with 
average column densities of 14.37 and 14.53, respectively. Finally, studies of extragalactic 
dwarf systems with strong current star formation show at least weak O VI absorption when 
the major exciting clusters are observed directly (see Heckman et al. 2001 for NGC 1705; 
Lebouteiller etal. 2004 for I Zw 36; the present work for NGC 625). A large sample of FUSE 
dwarf galaxy observations spanning a range of galaxy properties (star formation rate, outflow 
velocities, dust content, etc.) is needed to quantify the importance of O VI in the cooling 
portions of the outflows. 



4.2. Diffuse Molecular Gas 

Diffuse H2 is a ubiquitous component of the ISM of the Galaxy. Using Copernicus data. 
Savage etal. (1977) demonstrated a strong correlation between H I column or reddening and 
the presence of diffuse H2, with molecular gas usually found on sightlines with sufficient H I 
to allow self-shielding or sufficient dust to allow the efficient formation of H2. However, these 
early observations were limited to local sightlines due to instrumental limitations. The higher 
sensitivity of ORFEUS allowed H2 to be studied in the Magellanic Clouds (see, e.g., de Boer 
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etal. 1998; Richter etal. 1998). Most recently, FUSEhets greatly expanded the amount of 
observational data available on diffuse H2. However, even though FUSE oHers unprecedented 
spectral resolution and sensitivity (typical observations probe H2 column densities ~ 10^^ 
cm~^), there have remained few detections of H2 in extragalactic environments. Within the 
Local Group, diffuse H2 has only been detected in the Milky Way (see, e.g., Shull etal. 2000 
and numerous other studies), the Magellanic Stream (Richter etal. 2001; Sembach etal. 
2001), the Magellanic Clouds (Tumhnson etal. 2002), and M33 (Bluhm etal. 2003). As 
mentioned in § 1, extragalactic detections of diffuse H2 beyond the Local Group are even 
more rare. 

Dwarf starburst galaxies such as NGG 625 offer a unique opportunity to study diffuse H2 
in the ISM, since they are typically UV-bright, are forming stars rapidly (implying a sizable 
molecular reservoir), and usually present relatively low extinctions due to their low metal 
abundances. Even with these apparent observational advantages, however, stringent column 
density limits of log(N(H2)) < 15 in the low- J levels arc derived for the metal-deficient star- 
burst galaxies IZwl8 (Vidal-Madjar etal. 2000), NGC1705 (Heckman etal. 2001), Mrk59 
(Thuan etal. 2002) and lZw36 (Lebouteiller etal. 2004). The new sample of Hoopcs etal. 
(2004) also searches for diffuse H2 gas in NGG 3310, NGG 4214, M83, and NGG 5253; only 
M 83 and NGG 5253 show detectable H2 gas in absorption against the background UV con- 
tinuum of the massive stellar populations. 

Hoopes et al. (2004) interpret the low molecular fractions derived in their sample (typ- 
ically ^ 10~^; i.e., similar to the fraction found in the present study for NGG 625) to 
the stronger UV radiation fields of these galaxies. This stronger UV background will raise 
the minimum values of extinction and H I column necessary for diffuse H2 to remain in the 
ISM without being destroyed. They also note that some of the systems in their sample have 
been detected in GO tracer hues and that the inferred molecular masses exceed the derived 
diffuse H2 masses by many orders of magnitude. 

The detection of H2 absorption in NGG 625 highlights the importance of ISM geometry 
in the interpretation of these types of observations. As shown by Gannon etal. (2002), the 
dust distribution (and hence potentially productive areas of H2 formation) does not always 
follow the distribution of UV light. This implies that UV absorption experiments for diffuse 
H2 may not sample the sightlines expected to show the highest columns of molecular gas. In 
the case of NGG 625, we have a combination of sightlines into a moderately dusty starburst 
region (extinctions toward the massive clusters of up to ~ 1 magnitude in the V-band; 
see Gannon et al. 2003) that has a moderate UV radiation field (since no higher-level H2 
transitions are detected). Here, the potential for detection of H2 may be optimized, given 
the correlation between dust and H2 (see, e.g.. Savage etal. 1977). 
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The simplest interpretation of these H2 observations of nearby dwarf starburst and 
low-metaUicity galaxies appears to be one where the covering factor of diffuse H2 clouds 
is low. This implies that most of the molecular material is confined to relatively small, 
dense molecular clouds, as predicted in models of the low-metallicity ISM (see, e.g., Maloney 
& Black 1988; Norman & Spaans 1997; Pak etal. 1998; Bolatto etal. 1999). Since the 
strength of the UV radiation field appears to be an important factor both theoretically and 
observationally, it may be expected that sighthnes to the UV-brightest clusters may not 
show diffuse H2 in absorption. The complicated geometry of dust, neutral and molecular 
gas, and the current stellar populations will require a much larger sample to fully understand 
the nature of diffuse H2 absorption in extragalactic environments. 



4.3. Abundances in the Neutral and Nebuleir Gas 

In these data we find an apparent abundance difference between the neutral and nebular 
gas regions in NGC 625. Interestingly, the N/0 ratio is, within errors, identical in the neutral 
and nebular gas phases. We again emphasize that our O I A 1039.230 A line may be saturated 
(hence reducing the offset between N and O abundances), and that our absolute abundances 
are comparatively uncertain due to the large error on our measured H I column within 
NGC 625. However, taken at face value, the offset is large enough (~ 0.9 dex) to warrant a 
more thorough discussion of these values in other systems and the potential implications for 
the evolution of low-mass galaxies. 

As shown in Table 3, it appears that such abundance offsets between nebular and neu- 
tral gas are common in strongly star-forming dwarf galaxies studied to date with FUSE. 
Abundance differences of > 0.5 dcx have been found in NGC 1705 (Heckman etal. 2001), 
IZwl8 (Aloisi etal. 2003; but see also Lecavelier des Etangs etal. 2004 for an alternative 
treatment), IZw36 (Lebouteiller etal. 2004), and Mrk59 (Thuan etal. 2002). These abun- 
dance differences may be caused by FUSE sighthnes sampling two different components of 
the ISM of dwarf galaxies: lower-abundance halo gas, and higher-abundance gas nearer to 
the active star formation regions. 

A simple test of this "disk vs. halo" scenario can be performed by comparing the size 
of the abundance difference with the difference between the H I gas column seen in emission 
and the H I column that can only be foreground to the UV- luminous regions (i.e., the H I 

gas probed in FUSE observations). Our reasoning is simple: if the absorption arises in the 
disk of the galaxy, the column density should be comparable to the disk column density, but 
if the absorption arises in a halo, the column density should be smaller. We show in Figure 7 
a plot of the size of the offset between neutral and nebular gas (A([0/H]) = [0/H]neutrai ~ 
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[0/H] nebular) and the difference in H I columns as seen in absorption and in emission toward 

the same region (A(H l) = log(N(H i))fuse — log(N(H i))2icm)- Depending on the value 
chosen for IZwl8, there could be a trend of increasing A([0/H]) with increasing A(H l). It 
appears that each system under study (with the possible exception of lZwl8) shows some 
value of abundance offset between the nebular and neutral gas (see notes to Table 3 for 
more details). Each of these investigations is susceptible to various sources of error that will 
require a larger sample to overcome; however, if these results are supported by further data, 
it appears that many dwarf galaxies may have a halo of lower-metallicity gas that surrounds 
the actively star-forming regions as probed by emission line spectroscopy. 

This result has important implications for the chemical evolution of star-forming dwarf 
galaxies, since it implies that widespread enrichment episodes have preceded the current 
bursts that dominate the bolometric luminosity of these systems. Previous works have 
postulated a prompt "self-enrichment" scenario in the nebular gas region that would be 
evidenced by localized enrichment near current H ll regions (e.g., Kunth & Sargent 1986). 
This effect has been shown to be small, since local abundance variations are not seen in 
the nebular regions surrounding massive starbursts in low-mass galaxies (see Kobulnicky & 
Skillman 1996, Kobulnicky & Skillman 1997, Kobulnicky etal. 1997, Kobulnicky & Skillman 
1998, Legrand etal. 2000, and references therein). 

The present offset between neutral and nebular gas abundances would require a differ- 
ent enrichment scenario than "localized enrichment", since the absorbing columns are, by 
definition, foreground to the starburst regions. A potential scenario that would be consistent 
with these data is one where low-level star formation persists in dwarf galaxies for extended 
periods of time, allowing the oxygen and nitrogen abundances to be elevated in the ISM and 
efficiently mixed with the surrounding neutral gas. This low-level star formation rate would 
elevate the neutral gas abundances above the primordial value, and the higher nebular values 
could then be produced by more recent star formation in the current major star formation 
regions of these systems. The close agreement in N/0 between the nebular and neutral gas 
also points toward central creation and geometrical dilution over a long time period; other- 
wise the different production timescales for N and O become problematic. A larger FUSE 
sample of actively star forming dwarf galaxies would be most beneficial in addressing the 
abundance offsets between nebular and neutral gas phases. 

If the absorption spectra of dwarf galaxies obtained by FUSE are indeed sampling a 
halo of neutral gas and are not primarily produced in the disks, then they provide a very 
important probe of a virtually unstudied component of the dwarf galaxy ISM. Kennicutt & 
Skillman (2001) have emphasized the importance of measuring abundances in the neutral 
gas in the outer parts of dwarf galaxies. In spiral galaxies, it is well known that there are 
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chemical abundance gradients in the sense of lower abundances in the outer parts of the 
systems. However, in dwarf galaxies it is generally assumed that the entire H I disk has the 
same metallicity as measured in the H ll regions, and this is a very uncertain assumption. 
The physical basis for this assumption is the general uniformity of H II region abundances in 
dwarf galaxies (Kobulnicky & Skillman 1997 and references therein), and the inference that 
the whole H I disk is kept at a rather uniform chemical abundance by the rapid transportation 
of the metals in a hot phase of the ISM (Clayton & Pantelaki 1993; Tenorio-Tagle 1996). 
However, the H ll regions only sample the inner parts of the H I disk. In some dwarf 
galaxies, as much as 90% of the neutral hydrogen hes outside of the Holmberg radius (e.g., 
DD0154; Carignan & Freeman 1988; Carignan & Purton 1998). If dwarf galaxies do have 
chemical abundance gradients, then assuming that the chemical abundances arc constant 
overestimates the total metal content of the galaxy, and leads to a misinterpretation of their 
evolutionary status (e.g., artificially inflating the calculated effective yield). The edge-on 
orientation of NGC 625 implies that at least some of the absorption is occurring in the outer 
parts of the galaxy and thus providing a probe of this relatively unexplored ISM component. 

The extended mission of FUSE offers an ideal opportunity to test for this important 
evolutionary scenario. One would ideally seek a sample of luminous, metal-poor (< 10% 
Zq) dwarf systems that have low intrinsic H I foreground columns, low foreground and 
internal extinctions, along high-visibility sightlines not contaminated by intermediate- or 
high- velocity clouds. Deep integrations on such targets will provide sufficient S/N to allow 
inter-comparison of the columns derived from oxygen lines with different oscillator strengths 
throughout the FUSE spectral region. With the effects of line saturation ehminated, such a 
sample would allow the exploration of this potentially important ISM phase in dwarf galaxies. 



5. Conclusions 

We have presented new FUSE spectroscopy of the dwarf starburst galaxy NGC 625. 
These data allow a detailed investigation of multiple phases of the ISM of the galaxy, includ- 
ing the molecular, neutral, ionized and coronal gas contents. We use these data to study the 
kinematics of the ISM, the diffuse H2 content, and the abundances in the neutral gas phase. 

Our first major result is that O VI absorption has been detected in a relatively low- 
velocity outflow. This detection adds to the sample of extragalactic systems showing out- 
fiowing coronal gas. With these data alone we cannot constrain the efficiency of cooling in 
the coronal gas; however, comparing to the literature and given the strength of the O VI 
absorption feature, we suggest that radiative cooling is likely not a dominant mechanism in 
the loss of energy from the NGC 625 outflow. 
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Our second major result is that we have detected low rotational level transitions from 
the Lyman and Werner bands of diffuse H2 gas in NGC625. This is one of only a few 
dwarf galaxies outside of the Local Group to show diffuse H2 in FUSE observations. It 
is likely that the geometry of the stellar populations, dust, neutral and dense molecular 
gas all play a role in determining if H2 is detectable in such systems. We interpret the low 
molecular fraction, which is similar to that found in other dwarf starburst and low-metallicity 
galaxies from FUSE studies, to the low covering factor of H2 clouds, which are expected to be 
comparatively small and dense in these actively star-forming, metal-deficient environments. 

Our final important result is that we have found a potential abundance offset between 
the nebular and neutral gas phases of NGC625. While the absolute oxygen abundance is 
hindered by potential saturation and the intrinsic H I column is only constrained to ±0.3 
dex, the magnitude of the offset (~ 0.9 dex) suggests that it is hkely real. Interestingly, 
the N/0 ratio is found to be the same in both the nebular and neutral phases. Similar 
results have been found for other dwarf galaxies studied with FUSE, suggesting that dwarf 
galaxies may in general host a lower-metallicity halo of neutral gas that can only be probed 
by absorption line spectroscopy. A larger sample of systems will help to shed light on this 
important evolutionary scenario. 
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Fig. 1. — Position of the LiFl detector segment of the FUSE MDKS A" x 20" aperture over- 
laid on an HST/WY'PC2 V-band image (a) and on an HST/WY'PC2 continuum-subtracted 
Ha image (b). The field of view is ~ 875 x 471 pc at the distance of 3.89 Mpc, with north 
up and east to the left. The average position angle during these observations was 112° east 
of north, as shown. We estimate the positional accuracy of the sht to be ± 1", based on 
the nominal Ff/S'E pointing accuracy (~ 0.5" rms), the iK'r/WFPC2 pointing accuracy (~ 
0.5" rms) and the comparison to ground-based astrometric solutions. This slit placement is 
sensitive to the bulk of the UV flux from NGC 625. 
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Fig. 2. — Overview of our FUSE spectra, with the LiFlA spectrum shown in (a) and the 
LiFlB spectrum shown in (b). These data have been binned by 8 pixels (~ 0.05 A) for 
clarity. Note the numerous absorption lines detected throughout this spectral region. 
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Fig. 3. — Closer view of the spectral region surrounding Ly/5 and the O VI absorption 
features, from the LiFlA detector segment. Prominent absorption lines are labeled with 
both the Milky Way and NGC625 (systemic) velocities. The O VI absorption is strong 
(equivalent width ~ 0.2 A) and the velocity centroid is blueshifted with respect to the 
molecular and neutral gas absorption centroids. See further discussion in § 3.1, and also 
Figures 4 and 5. 
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Fig. 4. — Normalized line profiles of identified H2 lines arising within the ISM of NGC 625. 
Vertical hnes show Milky Way and NGC 625 velocities; these H2 lines are discussed further 
in §§ 3.1 and 3.2. 
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Fig. 5. — Normalized line profiles of identified neutral and ionized atomic absorption lines 
arising within the ISM of NGC 625. Vertical lines show Milky Way and NGC 625 velocities; 
these absorption lines are discussed further in §§ 3.1 and 3.2. 
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Fig. 6. — Profiles of Galactic and NGC625 H I Ly/9 absorption, (a) and (c) show the "low" 
and "high" continuum normalizations, and (b) and (d) show the resulting H I absorption 
profiles. Overlaid on (b) and (d) are profiles for the Galactic and NGC625 H I Ly/3 absorp- 
tion. For the Milky Way profile, three lines are shown at log(N(H l)) = 20.3 (solid line) and 
log(N(H l)) = 20.1, 20.5 (dotted lines). For NGC625, five lines are shown at log(N(H l)) = 
20.6 (solid), log(N(H l)) = 20.4, 20.8 (dotted lines), and log(N(H l)) = 20.2, 21.0 (dot-dash 
lines). From these profiles we deduce an absorbing H I column of log(N(H l)) = 20.5 ±0.3 
for NGC625, with the continuum normalization dominating the error budget. See § 3.3 for 
further discussion. 
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Fig. 7. — Differences between neutral and nebular gas abundances plotted versus the differ- 
ence between H I absorption column density as probed by observations and H I 21 cm 
emission column density, for 4 different dwarf starburst galaxies. In the key, A03 corresponds 

to the study of Aloisi ct al. (2003), and L04 corresponds to the study of Lecavelier des Etangs 
etal. (2004). Errorbars on A(H l) are difficult to estimate, and are shown at the 0.3 dex 
level. We find no obvious correlation between these two values; rather, abundance offsets 
between the neutral and nebular gas appear to be common in the dwarf galaxies shown here. 
To the left of the dotted vertical line at A(H l) = —0.3, one expects to probe "halo" gas via 
absorption line spectroscopy; to the right of this line, one expects only to probe gas close to 
the background sources within the disk. 
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Table 1. Basic Parameters of NGC 625 



Property 


Value 


Reference 


Distance (Mpc) 


3.89 ±0.22 


Cannon et al. (2003) 


Ms 


-16.3 


Marlowe et al. (1997) 


b, 


273.7°, -73.1° 




E(B-V) 


0.016 


Schlegel etal. (1998) 


12 + log(0/H) 


8.14 ±0.02 


Skillman et al. (2003b) 


Current SFR (Mq yr^i) 


0.05 


Skillman et al. (2003a) 


H I Mass (10** M0) 


1.1±0.1 


Cannon et al. (2004) 


^Helio (kms~^) 


413 ±5 


Cannon et al. (2004) 



^Galactic longitude and latitude. 
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Table 2. Important Line Parameters 



Line 


^0" 


-^NGC 626 


log(f-A)^ 


MsTGC 625 


WNGC625'' 


Iog(N)'= 


I.D. 


(A) 


(A) 




(kms-1) 


(A) 




H2 L R(0)(2-0)d 


1077.140 


1078.60 ±0.02 


1.099 


406 ±6 


0.07 ±0.007 


14.92 ±0.19 


H2 L R(l)(4-0)<1 


1049.960 


1051.38 ±0.03 


1.213 


406 ±9 


0.12 ±0.01 


15.01 ±0.15 


H2 L R(2)(5-0)'* 


1038.689 


1040.07 ±0.02 


1.234 


399 ± 6 


0.04 ±0.004 


14.34 ±0.19 


H2 L P(l)(5-0)d 


1038.157 


1039.51 ±0.02 


0.954 


391 ± 6 


0.05 ±0.005 


14.82 ±0.18 


H2 W R(2)(0-0)<^ 


1009.024 


1010.36 ±0.02 


1.198 


397 ± 6 


0.05 ±0.005 


14.59 ±0.18 


H2 L R(l)(8-0)'' 


1002.449 


1003.80 ±0.03 


1.262 


404 ±9 


0.09 ±0.009 


14.87 ±0.15 


C 11* 


1037.018 


1038.30 ±0.02 


2.088 


371 ±6 






N I 


1134.980 


1136.43 ±0.02 


1.674 


383 ±6 


0.10 ±0.01 


14.63 ±0.12 


N I 


1134.415 


1135.88 ±0.02 


1.512 


387 ±6 


0.09 ±0.009 


< 14.80 ±0.12 


I 


1039.230 


1040.55 ±0.02 


0.974 


381 ± 6 


0.24 ±0.02 


> 15.80 ±0.12"= 


VI 


1031.926 


1033.20 ±0.02 


2.136 


370 ±6 


0.20 ±0.02 


14.32 ±0.08 


Si 11 


1020.699 


1022.03 ±0.03 


1.234 


391 ±9 


0.19 ±0.02 


15.37±0.15 


P 11 


1152.818 


1154.31 ±0.02 


2.451 


388 ±6 


0.09 ±0.009 


13.54±0.16 


S III 


1012.495 


1013.79 ±0.02 


1.647 


383 ±6 


0.38 ±0.04 


15.04 ±0.15 


Ar I 


1066.660 


1068.10 ±0.02 


1.857 


405 ±6 


0.06 ±0.006 


14.20 ±0.20 


Ar I 


1048.220 


1049.58 ±0.02 


2.440 


389 ± 6 


0.12±0.01 


13.93 ±0.10 


Fe II 


1144.938 


1146.45 ±0.02 


1.978 


396 ± 6 


0.28 ±0.03 


14.72 ±0.16 


Fe II 


1063.176 


1064.61 ±0.02 


1.765 


404 ±6 


0.17 ±0.02 


14.73 ±0.15 



^AU atomic data arc taken from Morton (2003). All molecular data is taken from the H200IS website; 
see http://www.pha.jhu.edu/ stephan/h2ools2.html and McCandliss (2003). 

''Representative errors on the equivalent width values are ± 10%. 

"^Column densities derived from the apparent optical depth method; see Savage & Sembach (1991). 
'^The labels "L" and "W" denote Lyman or Werner bands, respectively. 

"The O I A 1039.230 A line may be saturated; without a second oxygen line for comparison, there is no 
empirical way to test the saturation of this line. We note this uncertainty and discuss it further in § 3.3.3. 



Table 3. Neutral vs. Nebular Abundances in Low-Metallicity Galaxies^ 



Galaxy 


Nebular 


Neutral 


Nebular 


Neutral 


Nebular 


Neutral 


A{H I) 


Reference(s) 




[N/H]'^ 


[N/H]- 


[0/H]b 


[o/ur 


log(N/0)'^ 


log(N/0)<= 


(dex)d 




NGC 625 


-1.07±0.13 


-1.80 ±0.34 


-0.47 ±0.07 


> -1.36 ±0.33 


-1.33 ±0.02 


< -1.17±0.15 


< -0.5 


1, 2, 3 


NOG 1705 


-1.47±0.14 


-2.16±0.24 


-0.45 ±0.07 


-1.23 ±0.22 


-1.75 ±0.06 


-1.66±0.11 


- 1.1 


4, 5, 6 


I Zw 18" 


-2.31 ±0.13 


-2.88 ±0.11 


-1.45 ±0.05 


-2.06 ±0.28 


-1.59 ±0.05 


-1.54±0.26 


< ±0.2 


7, 8, 9 


I Zw 18' 


-2.31 ±0.13 


-3.00 ±0.18 


-1.45 ±0.05 


-1.36 ±0.80 


-1.59 ±0.05 


-2.401°;! 


< +0.0 


7, 10, 9 


IZw 36 


-1.65 ±0.11 


-2.81 ±0.38 


-0.89 ±0.05 


-1.76 ±0.508 


-1.49 ±0.01 


-1.80 ±0.608 


-0.2'' 


11, 12, 13 


Mrk59 


-1.46±0.11 


-2.08 ±0.59 


-0.67 ±0.05 


-1.66 ±0.30 


-1.52 ±0.01 


-1.15±0.5 


N/Ai 


11, 14 



Referenees. — 1 - Skillinan etal. 2003b; 2 - This work; 3 - Cannon et al. 2004a; 4 - Lee & Skillman 2004; 5 - Heckman et al. 2001; 6 - 
Meurer etal. 1998; 7 - Skillman & Kennicutt 1993; 8 - Aloisi etal. 2003; 9 - van Zee etal. 1998; 10 - Lecavelier des Etangs etal. 2004; 11 - 
Izotov & Tliuan 1999; 12 - Lebouteiller etal. 2004; 13 - Stil & Israel 2002; 14 - Thuan etal. 2002 

'^The abundances here use the solar nitrogen value from Holvircger (2001) and the solar oxygen abundance from Asplund ct al. (2004). 

''Derived from nebular spectroscopy of the most luminous H II regions in the galaxies, since these will be most heavily weighted in UV 
FUSE measurements. 

''Derived from FUSE spectroscopy. 

'^Defined as the difference between the H I absorption column density toward the heavily-weighted UV sightlines as probed by FUSE, 
and the H I 21 cm emission column density (for the area nearest to the FUSE aperture): A(H l) = log(N(H i))fusb — log(N(H i))2i cm. 

''Neutral gas values from analysis of FUSE data by Aloisi etal. (2003). 

^Neutral gas values from analysis of FUSE data by Lecavelier des Etangs etal. (2004). 

8Calculated using Phosphorus as a tracer of oxygen abundance; see discussion in Lebouteiller etal. (2004). 

''Only peak H I 21cm emission column density value is available; see Stil & Israel (2002). 
'Detailed H i emission observations of Mrk59 are not yet available in the literature. 



